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1. Introduction 

Indirect CP violation in K decays, quantified by |£|, lias been playing an important role in 
finding the location of the apex of the unitarity triangle in the p-Tj plane and in constraining new 
physics, especially the structure of flavor changing neutral current in it. Experimentally |£| has 
been determined precisely as l£|=(2.233±0.015)x 10^^^ |jl|]. Within the standard model, |£| can be 
expressed as 

|£| = (known factor) x Ba:(/i) x f{p,f}), (1.1) 
where f{p,fl) is a known function of the Wolfenstein parameters, p and fj, and 

„ , . (K^\dY^{\-Y5)sdY^{l-Ys)s\K'^) 

The purpose of this work is to determine the parameter Bk with high precision using lattice QCD 
to give a strong constraints on p and fj through eq. ( |1.1[ ). 

As seen from eq. ([L^, the AS = 2 four-quark operator has the form of {V — A) x (V — A). 
This makes the lattice calculation of Bk much simpler if we take the overlap fermion formalism 
because the overlap fermions respect the lattice variant of chiral symmetry exactly at a finite lattice 
spacing, and as a consequence the mixing with operators with other chiralities is prohibited. This 
simplification makes a precision lattice calculation possible. 

We perform the calculation on a 16^ x 32 lattice using the RG Iwasaki action at j8 = 2.30. To 
accelerate HMC, we have introduced extra Wilson quarks and ghosts At a price for the accel- 
eration, the topological charge is frozen during the HMC evolution. Because of this, configurations 
are generated at a fixed topological charge 2 = 0. Six sea quark masses are taken in the range of 
[0.015,0.100] in lattice unit, which roughly corresponds to [1/6 x m^^^^ ,mf^^^] in physical unit. Our 
lightest pion mass is about 290 MeV, and gives nij^L ~2.7. The lattice spacing l/a = 1.67(2)(2) 
GeV is determined by ro = 0.49 fm in the Q = sector. The physical spatial volume of our lattice 
is about (1.9 fm)^. In order to study the topological charge dependence, we have also generated 
configurations at 2 = — 2 and —4 at msea=0.050. We have accumulated 10,000 trajectories for 
2 = and 5,000 for Q = -2 and -4. 

Calculations are done every 20 trajectories at each mgea- Six valence quark masses take the 
same values as those of sea quarks. All degenerate and non-degenerate mesons have been cal- 
culated. We employ Coulomb gauge for the gauge fixing condition except for the calculation of 
non-perturbative renormalization constant, in which Landau gauge is used. Low-mode averaging 
is implemented for all correlation functions, which substantially improves statistical signals. 

2. Method and results 

Two-point functions are obtained in the standard way with a wall source at ?src and a point sink 
at t. We repeat this calculation four times with fsrc=0, 8, 16, 24, and take an average over them. 
The axial two-point function is defined by and fitted to 

cf^Tit) = UO\A,it,x)iAr'\0))UO) - ^2^fp'np [e''"^' + e'"^^'-''^) , (2.1) 
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where 



A4{t,x)=qi{t,x)y4Y5q2it,^), = \^qi{t,x)j 74/5 \^q2{t,y)j , (2.2) 

q'2ix) = [1 -Dov/(2mo)]^2W, Zjf = {P\ £^i(0,f) 7475^2(0,0) |0), (2.3) 



fptrip = (0|A4|f ), mo = 1.6 and ^3 = 16^'. nip and ZJ^^^fp are extracted by a correlated fit. 

In the calculation of three-point functions, the meson (anti-meson) interpolating operators with 
wall source are put at fixed time slice ti {12) while the position of four quark operator t is varied. 
Three-point functions are repeatedly calculated with (f2,?i)=(8,0), (0,24), (16,8), (24,16), (16,0), 
(24, 8). These six pairs of t\ and t2 are classified into two sets by the time separation, \t2 — 1\ |=8 (24) 
or 16, which we call set A and B, respectively. Within each set, all the three-point functions are 
equivalent after proper translation in the time direction, so they are averaged after shifting. As for 
set B (with \t2 — 1\ |=16), two equivalent regions, < ? < 16 and 16 < ? < 32 are further averaged. 

The three-point function is defined by eq. ( p!4| ) and fitted to the form in eq. (2.5): 



C^iy^h^tA) = I(0| {Ar\t2)y dlXit,x) {Ar\h)y |0) 



(2.4) 
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(2.5) 



'?i7ai(1 ~ 75)^2 ^i7ai(1 ~ 75)'72 AS=2 four-quark operator defined on the lattice, and 
Nt = 32. The first term in eq. ( |2.5[ ) contains the hadron matrix element relevant to the calculation of 
Bk. Since the time direction of our lattice is not so large, we add two additional terms to represent 
an excited state contamination and a contribution wrapping around the lattice. The mass of the 
excited state mpi appearing in the second term is extracted from the point-point pseudoscalar two- 
point function. We confirmed that nipi are consistent with the experimental value of 7r(1300) in 
the chiral limit within the error. The third term, expressing a wrapping contribution, contains a 
two-meson system, and the energy shift Ap = ftotai — 2mp is extracted from the fit. 



We simultaneously fit two sets of three-point functions to eq. (2.5) with mp fixed to the value 
extracted from the two-point function. As seen in Fig. |l|, ^-dependence of the three-point functions 



are well described by eq. ( |2.5| ). Then the lattice B-parameter Bp' is obtained by 
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Figure 1: Three-point functions. Data and fit results are shown. The fit range is [fmm,fmax] — [14,26] and 
[6, 10] for set A and B respectively. 
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Figure 2: Sea quark mass dependence of B^'^(2GeV). The different symbols denote different valence quark 
mass: 0.015-0.10 from bottom to top. Only the data consisting of degenerate valence quarks are shown. 
Lines are the results from a linear fit and just a guide to eyes. 



where the first and second factors are obtained from the two- and tiiree-point functions, respectively. 
The fit range dependence of Bp was studied, and found to be stable. 

We adopt the RI-MOM scheme to calculate the renormalization factor. Following the standard 
method, we obtain a preliminary result 

Z|^i = 1.217(6), z|f(2GeV) =0.862(4) (the error is statistical only). (2.7) 

3. Test of NLO ChPT and extraction of Bk 

We first test whether the quark mass dependence of Bp is consistent with the NLO partially 
quenched ChPT (PQChPT) prediction, or to which quark mass the prediction describes data well. 
In the test, we only use data points which satisfy wzgea < "^valence for the reason described below. 

Figure ^ shows the sea quark mass dependence of Bp, in which clear dependence is not seen 
except for the region with msea > 'Mvaience- In Ref- [01, the finite volume effects to Bjc were studied 
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test with NLO ChPT with f fixed 
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Figure 3: Test with the NLO ChPT formula. The different symbols denote the different sea quark mass. 
The lines correspond to those in msea = 0. In the figure legend, [0.015,0.100] denotes the range of the sea 
and valence quark masses used in the fit, for example. 

to NLO in the framework of PQChPT, and found to become more significant when wigga > wivaience- 
While the size effect found in Ref. ^ is tiny, it is pointed out in Ref. [Q] that the NLO estimate 
significantly underestimates for mj^ and fj^. For example, the NLO estimate of the size effect to /;;: 
gives about 2 % correction at our lightest unquenched point while the inclusion of NNLO gives 4-5 
%. Motivated by these observations, we include the data point in the fit only when mgea < 'Wvaience- 
The test is made using data consisting of degenerate quarks. Bp is fitted to the NLO PQChPT 
formula Igll, 



Bp 



Bj 



6mp , / m 
1 — -T-. — ^TTT In 



{bi -bj,)m\ + b2m^^, 



(3.1) 



where m^^, ~ B(){m^sa + "Jsea) and the free parameters are Bp, /, {b^ — b^) and bi- The fit results 
are shown in Fig. |3| (left). /, the tree level pion decay constant with / ~130 MeV normalization, 
and ;^^/dof are also shown for each fit range. While all fit ranges tested give acceptable x^/dof, 
/ monotonically decreases as the fit range is made narrower, /'s obtained from the two narrowest 
ranges are consistent with each other within one standard deviation, and its value is consistent with 
a naive expectation / = 100 ^ 130 MeV. Fixing / to 110 MeV [g], we obtain Fig. | (right). The 
;^^/dof values suggest that while the PQChPT formula does not apply for the two heaviest data 
points, the data for lUq < 0.050 (roughly corresponding to half strange mass) are inside the NLO 
ChPT regime. 

To extract Bk, we fit the data of both degenerate and non-degenerate quarks to the following 
formula Iglig], 



B 



12 



12 



{4nfy 



+ — 2 1""^12 



In 



m 



12 



+ 2 In 



+b\ + mj 1 




1 + C3 m j2 + C4 m j2 ' 



(3.2) 
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Figure 4: m^j dependence of Bp'^(2GeV). The different symbols correspond to the different m^, 
lines represent Bp extrapolated to msea = Wvi = "^uf/**- The vertical line denotes the position of physical niK- 



where m? ~ Bq ("^v; + ni^j) and m,,,- denotes a valence quark mass. The last two terms in eq. (3.2) 



are added to describe the data in the heavy region. The fit is performed with four data sets, each set 
including data from lightest three, four, five and six sea quarks. 

The fit results are shown in Fig. ^. The solid line is the one in which the lighter valence mass 
(wivi) and the sea quark mass (mgea) are extrapolated to the physical u^d mass {n^'^^^)- Interpolating 
to physical m^, we obtain B^^{2 GeV)=0.533-0.523 depending on the data sets used. As our 
preliminary result we take the result using four msea data, and obtain 



Bf^{2 GeV) =0.526(9) 



(3.3) 



where only the statistical error is shown. 

Since in the above fit we did not include the data which could have potentially significant finite 
size effect, the fit result is expected to be under control. As a conservative upper bound of the finite 
size effect, we take that of f^, and add a 5% error. 



4. The effect of fixing topology 

To estimate the effect of fixing the topological charge on Bp, according to the studies in 
Refs. [|8|, ^, we assume it to be 



l-7^h (4-1) 



(47r/)2 {Q2) V {Q 

where (Q^) = XiVa. ~ 10 at ntq = 0.05 [10|. This is motivated by an observation that the most 
significant -dependence of the physical quantities is that of pion mass, and other quantities are 
affected through it. Then, the correction to the Q = result is estimated to be 1.4% at mg=0.05, 
and the difference between 2 = and —2 (—4) to be 0.6% (2.2%). Since the size of the statistical 
error for Bp is about 2%, one does not expect to see clear Q dependence of Bp. In Fig. |5|, Bp 
at nisea = 0.05 from three Q are compared, where only the data of degenerate quarks are shown. 
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Figure 5: Compai-ison of Bp^(2GeV) at msea=0.05 with three different Q. 

We could not observe any systematic Q dependence which is statistically significant. Thus the 
assumption eq. seems to give a reasonable or even conservative estimate. We will quote 1 .4 
% as a crude estimate for the systematic error due to fixing topology. 

Since the calculation is made only at a single lattice spacing, it is difficult to estimate system- 
atic uncertainty due to scaling violation, though this is expected to be under control as no 0{a) 
error is present. The renormalization factor Z^^ may have a sizable systematic error as well. The 
study to estimate all these errors is in progress. 

Numerical simulations are performed on IBM System Blue Gene Solution at High Energy 
Accelerator Research Organization (KEK) under a support of its Large Scale Simulation Program 
(No. 07-16). This work is supported in part by the Grant-in-Aid of the Ministry of Education (No. 
17740171, 18034011, 18340075, 18740167, 18840045, 19540286, 19740160). 
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